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CONCENTRATION ALONG GEODESICS FOR A NONLINEAR STEKLOV PROBLEM 

ARISING IN CORROSION MODELLING 

CARLO D. PAGANI, DARIO PIEROTTI, ANGELA PISTOIA, AND GIUSI VAIRA 


Abstract. We consider the problem of finding pairs (A, u), with A > 0 and u a harmonic function in a three 
dimensional torus-like domain "D, satisfying the nonlinear boundary condition 9„u = A sinhu on 8T>. This type 
of boundary condition arises in corrosion modelling (Butler Volmer condition). We prove existence of solutions 
which concentrate along some geodesics of the boundary dT> as the parameter A goes to zero. 


1. Statement of the problem and main results 

In mathematical modeling of electrochemical corrosion a suitably defined galvanic potential satisfies an elliptic 
equation (namely, the Laplace or Poisson equation in simplest cases) in a given domain T >, whose boundary is 
partly electrochemically active and partly inert. In the inactive boundary region, the current density flow is 
of course zero, but in the active part it is modeled (by interpolating experimental data) by a difference of two 
exponentials according to the so-called Butler-Volmer formula (see [2] for a detailed discussion of the model). 
Then the resulting mathematical problem consists in solving an elliptic equation complemented with a boundary 
condition of Neumann type, namely 


( 1 . 1 ) 


d n u(y) = A y(y) [e au ^ - e"* 1 " 01 ^], y € dV 


Here n is the outward unit normal to dT>, u is the surface potential, a £ (0,1) is a constant depending on the 
constituents of the electrochemical system, y(y) is a non negative bounded function which distinguishes between 
the active and the inert boundary regions (typically fi(y) is the characteristic function of some subset of dV) 
and A is a positive parameter. 

Due to the exponential grow of the nonlinear boundary term, this problem has beeen studied (usually for 
Laplace or Poisson equations) in two dimensions: see [II] . [7], [8], [Sh 

In the physically relevant three dimensional case, little is known about existence and properties of solutions (see 

m- on])- 

A case which presents some interest for applications arises when the body T> is a three dimensional annular 
shaped domain, namely it can be represented in the form 


(1-2) V = | (yi,J/ 2 , y 3 ) G R 3 ( 

and H is a smooth bounded domain in R 2 such that 


y/vi+ylva) G 


(1.3) fl C {(xi,x 2 ) G R 2 | x\ > 0} 

Clearly, the domain T> is ^-invariant for the action of the group Q given by 


g(yi, 2 / 2 , 2 / 3 ) = (^( 2 /i, 2 / 2 ), 2 / 3 ) 


where g £ 0(2) (the group of linear isometries of R 2 ). 

The chosen geometry aims at modellizing the corrosion of torus-like bodies. Then, we consider the boundary 
value problem 


r Au = 0 in T> 

\ d n u = A sinh u on dV 


Note that we have chosen in ED a = 1/2, fi = 1, and we wrote u instead of u/2 for the harmonic potential. 
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In such situation it is natural to look for solutions which are (/-invariant, i.e. they are axially symmetric functions 
of the form 


(1-5) u(yi, 2 / 2 , 2 / 3 ) = u{x 1 , x 2 ) 

where X\ = ( y\ + y I) 1 / 2 and X2 = 2 / 3 - 

Now an easy computation shows that (El is equivalent to the following problem for u : 


( 1 . 6 ) 


div(xiVu) =0 in 
d v u = A sinhu on dLl 


Thus, we are led to study the more general anisotropic two-dimensional problem 


(1.7) 


f A a u = 0 in D 
( d v u = A sinh u on dLl 


with a G L°°(f2) such that 

(1.8) 0 < ao < a(x±,X 2 ) < ai < 00 a.e. in f2, 


A is a positive parameter and A Q zi is defined as 


(1.9) 


A a U 


, . div(a(a;)Vu) 
a(x) 


Au + V log a(x) • Vu 


for every u G U 1 (17). 

We remark that problem m corresponds to choosing a(x 1 , 2 : 2 ) = x\. 


The main goal of the present paper is to extend to the anisotropic problem El the principal results obtained 
for a = 1 both concerning existence of multiple solutions and the limiting behaviour of the solutions for A —> 0 + . 

In what follows we first (section 2) prove existence of global multiple solutions of problem (11.711 and, as a 
consequence, of problem El- This is done by a mild modification of the variational approach used in [9]. We 
prove the following result: 

Theorem 1.1. Let LI C R 2 be a Lipschitz domain and let a G L°°(Ll) such that (11.81) is satisfied. Then for 
every A > 0 there exist infinitely many solutions (in H l (Ll)) of problem (11.71) . 

Then, a simple Corollary states the existence of solutions of problem El: 

Corollary 1.2. For every A > 0 there are infinitely many axially symmetric solutions to problem El- 

Moreover, in subsection 12.31 a first approach to the description of the behaviour of the solutions for A —» 0 + is 
considered. For a = 1 this study was done in [7]; there the authors prove, for a bounded C 2, “ domain C R 2 , 
that all solutions have energies that are of order log(l/A); the limit of the boundary flux along subsequences 
is a sum of Dirac masses located at a finite set of points, potentially accompanied by a regular part of definite 
sign. These results are a generalization of those obtained in |3] by studying the behaviour of explicit solutions 
of the problem in the unit disk. We show that the same kind of results hold for the anisotropic case, and we 
sketch the main arguments of the proofs. 

In [6] , considering again the case a = 1, the authors prove that in any domain Q there are at least two distinct 
families of solutions which exhibit exactly the qualitative behaviour of the explicit solutions found in [3] , namely 
with limiting boundary flux given by an array of Dirac masses with weight 2n and alternate signs. We will show 
in section [3J that, for the anisotropic problem EH), the situation is substantially the same as depicted in ©. 
In that case the arguments used to extend the results of [6] to the anisotropic case are in some part different 
from those given in ;> 6j. 

In order to state our result let x = ( Xi,X2 ) G Ll, y = ( 2/1 , 2 / 2 ) G dLl and let G a (x,y) be the Green’s function for 
the Neumann problem 
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( 1 . 10 ) 


-(A a ) x G a (x, y) = 0 in fl 

dG a 0 r , 2 tt 
~Q t — = 2 7 r<>y(a;) - 7— Z 7 Z\ a yy> 


fan a ( x ) 


on 9S1 


[ G a {x,y) = 0 

Jan 


and define H a (x,y) as the regular part of G a {x,y ), namely 

(1-11) H a {x,y) = G a {x,y) - log ,- 

F - yr 

We say that £ G <9S2 is a C 1 —stable critical point of a|an if Vpna(£) = 0 and the local Brouwer degree 
deg (Vanci, B{f, p) n<9f2, 0) is well defined (for p small enough) and it is not zero. It is easy to see that any 
strict local minimum point or strict local maximum point or non degenerate critical point of a|an is a C 1 —stable 
critical point. 


Now, the main result is the following 


Theorem 1.3. Let Q. C ffi 2 be a smooth domain and a G C 1 (n). Let £i,£ 2 G dfl be two different C 1 — stable 
critical points of a\an- Then there is Ao > 0 such that, for 0 < A < Ao, there is a sign-changing solution u\ of 
the problem m of the form 


( 1 . 12 ) 


U\(x) = log 


2/h _ 

(£1 + Apii'i)| 2 



2/r 2 _ 

(£2 + A/i 2 ^ 2 )| 2 


+ 0 ( 1 ) 


where V\ and ^ 2 denote the unit outer normals to dLl at the point £1 and £ 2 respectively, and the parameters pi 
and /i 2 are explicitly given by 

Ml = I e ffa(Cl,ei)-G a (6A2) j = I e iJa(?2,6)-G a (e2,«l)_ 


In particular, the solution u\ concentrates positively and negatively at the points £1 and £ 2 respectively, as A 
goes to zero. 


According to the previous discussion, the corresponding result for problem (11.41) reads as follows: 

Corollary 1.4. Let fl C K 2 be a smooth domain and a G C 1 )!!). Let £i,£ 2 G dfl be two different C 1 — stable 
critical points of a(x ) = x\ restricted on dfl. Then there exists Ao > 0 such that for any A G (0, Ao), problem 
JO) has a sign-changing axially symmetric solution ua which concentrates positively and negatively along two 
geodesics of the boundary dT> which are the Q— orbits of f \ and £ 2 respectively, as A goes to zero. 


The proof of Theorem 11.31 relies on a very well-known finite dimensional procedure and it is carried out in 
Section [3] We shall omit many details of this proof because they can be found up to minor modifications in [Bj, 
where problem Gnu has been studied with a(x) = 1. We only compute what cannot be deduced from known 
results. It must be mentioned that in our case the anisotropic term a allows to solve the reduced problem by 
finding the location of the concentration points in terms of critical points of a itself. 

Our result suggests that it should be possible to find solutions which concentrate along two geodesics of the 
boundary of a more general torus-like 3-dimensional domain which is not necessarily axially symmetric. 


2. Existence and multiplicity of solutions 

The content of this section almost repeats the arguments developed in [9j : recall that we have chosen the weight 
function p{x) appearing in (11.11) identically 1, and that we deal with the operator A a instead of the Laplacean; 
therefore many statemets are simple rephrasements of the corresponding statements given in [9] and are not 
given here in full details. Let us first recall that the approach to the non linear problem relies on the solution of 
a related linear Steklov eigenvalue problem on the boundary. We summarize here (without proofs) the crucial 
results about this problem. 
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Let n C 

H^ny. 

( 2 . 1 ) 


be a bounded Lipschitz domain and consider the following linear Steklov eigenvalue problem in 

J div(a(a;)Vw) =0 in fi 
I d„u = Ait on dLl 


It is easily seen that for A > 0 the solutions to (EH> belong to the subspace H\ C H l (Vl) defined as follows: 


( 2 . 2 ) 


Hi = ju e H\n), J au = o|. 


It can be shown (by a classical reductio ad absurdum argument, see e.g. 0) that in H\ the Dirichlet norm 
f n |Vu| 2 is equivalent to the H 1 norm and that (12.11) is equivalent to the following variational problem: 

Find u £ H^, u 0, such that 

(2.3) / aV«V^ = A / amp 

Jn Jan 

holds for every p £ H 1 . Furthermore, the expression 


(2.4) 


lull" = / a |Vw 
In 


an 


defines an equivalent norm in U 1 (12). We will consider the scalar product in IF 1 (r2) associated to this equivalent 
norm; by a slight adaptation of the proof in [5], theorem 2.1, we have the following result : 

Proposition 2.1. Problem (12.31) has infinitely many eigenvalues 

0 < Ai < A 2 < ... 

each of finite multiplicity and such that |A„| —>• + 00 . The eigenvalue Ao = 0 corresponds to the constant 
solutions of the homogeneous Neumann problem. Moreover, we can take all the eigenfunctions v n , n > 0 
orthogonal and normalized with respect to the scalar product associated to the equivalent norm m and the 
following decomposition holds: 

(2.5) H 1 = Hi © V a 

where the subspace V a is spanned by the eigenfunctions v n (note that V a contains the subspace of the constant 
functions). 

Remark 2.2. Global regularity of the eigenfunctions depends on the weight a(x) and on the regularity of the 
boundary dfl. It can be shown that on a Lipschitz domain If the eigenfunctions are bounded and continuous on 

O. 

2.1. Main estimates. Let us consider now the non linear problem (El); we define the even functional 

(2.6) E\(u) = i f a|Vw| 2 — xf a(coshw —l) 

2 Jn Jan 

where u £ R 1 (n). Note that 

(2.7) 

where, for every q > 1, 


cosh it — 1 = — + u^hfu) 


h : H^fl) ->• L q (dfl) 

is bounded (see [7], lemma 2.1). 

Define further 

(2.8) S R = {u£H 1 : ||n|| a = i?}; 

where || || a is the equivalent norm defined in (12.41) . Then, we have 

Lemma 2.3. For every A > 0, there exists R > 0 and a closed subspace V + C H 1 ^) with codimV + < 00 such 
that 

E\{u) > c 0 > 0, 


for every u £ Sr D V + . 
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Proof. By (12.71) we have 

(2.9) E\(u) > ^ [ a |Vm| 2 - ^ [ au 2 - A [ a\u\ 4 h(u), 

z Jn 1 Jon Jon 

and for every q > 1, p = q/(q — 1), the integral in the last term can be bounded as follows 

, i/p 


( 2 . 10 ) 


I dQ 


a \u\ 4 h(u) < ||a|| L =o (an) ||/i(it)|| L , ( an)( [ M 4p ) < C 

y Ja n J 


<C\\u\\l = CR\ 


^k — ^ ^ ^ k-\-r 

V + = H l 0 © span„> fe+r { v n }, 


Let us now consider the quadratic part of the functional. Since the sequence of positive eigenvalues \ n of the 
linear problem is unbounded, there exist non negative integers k, r, such that 

(2.11) 

Then, we set 

(2.12) 

For every u £ Sr D V + we get 

(2.13) l[a\Vu\ 2 -^f au 2 > - y —) j a|Vu| 2 = ^(l - 

z ./Si ^ Jd n z v i \ Afe+ r / 

Then, the lemma follows by taking R small enough. □ 

We are now going to construct closed, finite dimensional subspaces, V~ C H 1 ^) such that: 

• dim V~ > codim V + \ 

• E\{u) < Coo < oo for every u £ V~ 

To this aim, let v ni , 1 < i < l be any finite sequence of l eigenfunctions, with l > codim V + , corresponding to 
the eigenvalues 

(2.14) 

Let us define: 


^ni — ^n,2 — — ^ni • 


(2.15) V =span 1 < i < i {c„ i }, 

The next lemma provides the key estimates at infinity on the functional (12.61) . 

Lemma 2.4. Let V~ be defined by (12.151) . Then, for every A > 0 we have E\{u) < 0 for any u £ V~ with 
large enough norm. As a consequence, there exists c aa < oo such that 


E\{u) < c 0 


Vit £ V~ 


Proof. We first assume A ni > 0. For notational simplicity, from now on we set v Ui = Vi, A ni = Ai > 0 
(1 < * < l). Thus, we can write any u £ V~ in the form 

i 

U = ^ ( tiVi 
z=l 

Recalling remark [?~21 u is a bounded continuous function on 0, so that sinh u £ H 1 ( fl) and by proposition 12.II 
the variational equation 


a Vcj V sinh u = A j 


i an 


avj sinhu 


(2.16) 

holds for every j, 1 < j < l. Multiplying by tj and summing up from j = 1 to l, we find 

(2.17) / a |Vu| 2 cosh(?/) = NJ Ajtj / avj sinhu 

JQ j=l J d£l 

On the finite dimensional space V— the functional (12.61) takes the form 

E\{u) = f(ti, ... ,ti ) = ^ f ~4r I a ^ v i ' - A f a (cosh u — l) 


i,j =1 


an 
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(2.18) a(coshu-l), 

2 i= i 

where we used orthogonality and normalization of v% with respect to the inner product defined by the equivalent 
norm (HU). Then, we have 

i i 

A jtj — A A jtj 

3 =i 3=1 


i an 


avj sinhu 


( y~! 


3 = 1 


(2.19) 


< A; y^ t 2 — X a |Vu| 2 coshu = / a |Vu|' 
J=1 J n J n 


A; — A cosh u 


Now, it can be proved (see [8], lemma 3.5) that the last term is strictly negative for ||u|| a = y/tf + ...t'f large 
enough. But the first term is the derivative of the function / along the curves 


h = cie Als , = cie x,s 


s G 


(orthogonal to the hypersurfaces Ai t\ + ... + A it 2 = constant); hence, for large yjt\ + ...t'f, the function / is 
strictly decreasing along these curves. 

We conclude that f(u) < 0 for u G V~ with ||it|| Q large enough; since / is continuous and V~ has finite 
dimension, we have 

sup f(u) = Cqo < °° 
u&v- 

We are left to show that we may allow A rai = 0 in (12.141) . This can be proved by the same arguments as in the 
proof of lemma 3.5 of [9]. □ 

2.2. Proof of theorem 11.11 To prove existence and even multiplicity of solutions to problem dE3 for every 
A > 0, we will apply the Symmetric Mountain Pass Lemma (see m thm. 6.3); thus, we need to show that the 
functional (HU) satisfies the Palais-Smale condition; to this aim the following estimate plays a key role: 

(2.20) 4(coshu — 1) < usinhu + u 2 

Proposition 2.5. Let z m G H 1 (fl) be a sequence such that E\(z m ) —> c and E' x (z m ) —> 0 in H 1 { ft)'. Then, 
the sequence z m is bounded and the functional m satisfies the Palais-Smale condition. 

Proof. Assume by contradiction (considering a subsequence if necessary) that ||z m ||o —> +oo and define 
tm = ||~m||a, Um = t" 1 Zm- Substituting in the condition E' x (z m )u = o(l)||u|| a , we get 


( 2 . 21 ) 


iVUmVu- A 


f sinh(t m u m ) 

/ a - - - U = 0[\)\\u\\a/t r 

I dQ tm 


Since u m is bounded in H 1 ( 12), there is a subsequence (still denoted by u m ) such that u rn converges weakly in 
H\n) and u m |an converges strongly in L 2 (<9fl); it can be proved (see [8], proposition 4.2) that u m |an —>• 0 a.e. 
By choosing u = u rn in (12.211) we get 


-A 


(2.22) / a |Vr 

J n 

On the other hand, from E\(z m ) ->cwe also get 

(2.23) 


f sinh(t m u m ) , , . 

/ a ---w m = o(l/t m ), 

/ an tm 


\ ( a |Vu m | 2 - A / 
z J n J as 


c o s h (1 rn u rn ) 1 


t an t m 

By comparison of (12.221) and (12.231) and by taking account of (12.201) we find 


= o(i/t 2 m ) 


0 = A / a 
Jan 


cosh(t„ 


- 1 


sinh(t m u m ) 


tl 


tr. 


~Um 


f cosh(t m u m ) - 1 

< —2A / a--b A 

J 80. 


an 


+ o(l/t m ) < 
o(l/t m ) 
























CONCENTRATION ALONG GEODESICS FOR A NONLINEAR STEKLOV PROBLEM 


7 


By recalling that u m —> 0 in L 2 (d£i) we finally get 


(2.24) 


0 < -2A / a 

Jan 


COSh (t m U m ) - 1 


+ o(l) 


By the above relation and again by (12.2311 we conclude 


a |Vu m | 2 -A 0 


so that ||u m || 0 —> 0, thus contradicting ||u m || a = 1- Then, the norm sequence ||z m || a is bounded and the same 
holds for \\z m ||. We can write 

Zm — C m 4“ Z m 

where c m is a bounded sequence and z m is bounded in H\ (12) (see definition (12.21) ). 

Now, the linear map L : —> ff 1 (12)' 

L(u)ip = / aVitVt^ 

Jn 

is boundedly invertible (by the Lax-Milgram theorem) while the operator 


T(u)ip = / a sinhuty? 


an 


maps bounded sets in iJ 1 (12) to relatively compact sets in iJ 1 (12) / (the result follows by an obvious extension 
of the arguments in [Tj , Lemma 2.1 and Lemma 2.2). By standard results [12] Proposition 2.2, it follows that 
z m is relatively compact in then, by the above decomposition, the same holds for z m in iJ 1 (12). □ 


Proof of Theorem ED' By lemmas \2 .31 and [2~T1 for any positive integer m there exist two closed subspaces V + , 
V~ of H 1 ^) with dimP - — codimP + = m, and positive constants R , Co, c Q0 (the last one depending on to) 
such that: 


a) E\{u) > co Vw £ W, |M|a = R', b) E\(u) < Mu G V . 

Then, by theorem 2.4 of [2], the functional E\ possesses at least m distinct pairs of critical points, corresponding 
to critical levels Cfc(A), k = 1,2, ...,to, given by 

(2.25) Cfc(A) = inf sup-E;,^) 

ueA 

where is the set of closed, symmetric sets A C such that 7 *(A) > k and 7 * is the pseudo-index related 

to the Krasnoselski genus 7 and to the subset Sr fl V + (see [5], definition 2.8). 

Moreover, 

Co + Cl 7 C2 + • ■ • + Cm + Coo ■ 

Since this conclusion holds for arbitrary to, we get infinitely many critical points; hence, problem m has 
infinitely many solutions in H 1 { fl). By standard regularity results [4], if f l is smooth and a £ C°°(fl), we have 
u e 

Remark 2.6. In the degenerate case Ck = ■■■ = Ck+ r = c (with k > 1 and k + r < m) it was shown in [2] that 
7 (K c ) > r + 1 > 2 , where I\ c is the set of critical points at level c; since a finite set (not containing the origin) 
has genus 1, it follows that E\ has infinitely many critical points at level c. 

Remark 2.7. The results of existence and multiplicity obtained in [7] correspond to the case a = 1 of Theorem 
11.11 above. 

From the discussion given in the introduction, the existence of solutions of problem (ED stated in Corollary 
11.21 easily follows. 
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2.3. Estimates and limits for the variational solutions. As discussed in the introduction, we now prove 
lower and upper bounds for the variational solutions obtained in the previous subsection; more precisely, we 
show that a branch of solutions corresponding to any of the critical levels c&(A) blows up in energy (as well 
as in the Dirichlet seminorm) at the rate log(l/A) for A —>• 0 + , while the corresponding normal currents stay 
bounded in L 1 {dCl). 

We first establish the lower bounds; by observing that a solution u to ed satisfies 


0 = / ad„u = \ / a sinh u 
J dQ J dQ 

we may write u = u° + s, s £ K, where u° solves the problem 

( A a u° = 0 in Cl 


(2.26) 

and 


d v u° = A sinh(u° + s ) on d fl 


/ au° = 0 , / a sinh(u° + s) = 0 

J dQ J dQ 


By elementary calculations, the second identity yields 
(2.27) 


s = s ( u °) = ilog%^_ 
2 


Now, following [7i and by recalling the bound (11.81) . one first proves the inequality 

|s(u°)| <C'i+C 2 ||Vu 0 || 2 


(with positive constants C\, C 2 depending only on Cl and on a) and subsequenly the following lower estimates 
for the solutions to ED 

(2.28) E x (u)> Alog(i)-B, ||Vu || 2 > Alog(i) - B 

where E\ is the energy functional defined in (12.61) and the positive constants A, B are independent of A and u. 
By considering the variational solutions given by theorem 11.11 we now provide upper estimates of the critical 
values Cfc(A) (12.251) . We first remark that any finite dimensional subspace V with dirnE > codimV' + ( V + being 
the subspace defined in ( 12 . 121 ) ) satisfies V € £*, for k < dimE — codimE + (this is related to the so-called 
intersection lemma, see 0 . lemma 6.4 and the proof of theorem 2.4 of 0 )- Moreover, since for A —» 0 + we can 
take V + = H\, that is the subspace of codimension 1 defined in ( 12 . 21 ) . we may assume dim V = k + 1. 

Now, let us choose k + 1 distinct points on the boundary dQ and, for e > 0, define 

V e = span {— log(e 2 + \x- aj | 2 ), j = 1,2,..., k+1} 

By the estimates of lemmas 3.4 and 3.5 in [7; (with obvious modifications due to the expression (12.61) of E\) 
one finds, for A sufficiently small and by taking e of order A : 

max£^(u) < CAlog^) 
uey A 

for some constant C* depending only on k, a and 0. Since by the previous remark we have V € G £&, the bound 

(2.29) Cfc(A) < C* log(^) 
follows. 

Finally, the estimate (12.291) together with an elementary lemma in integration theory (lemma 3.6 of [?]) yield 
the above mentioned bound on the normal current: 


Proposition 2.8. Let A > 0 be small enough so that (12.291) holds. For k > 1, let Uk,x be the variational 
solutions obtained in theorem o with E\(uk) = Cfc(A). Then , there exists a constant D * depending only on k, 
a and Cl, such that 


f 

dUk, A 

Jan 

dv 



sinh(u fciA )| < D* 


(2.30) 
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The proof follows by a trivial modification of the proof of Corollary 3.7 in [7]. 

By exploiting the previous estimates it is possible to describe the behaviour of a sequence of solutions u\ n = 
Uk, a„, where k > 1 is fixed and X n —> 0. The results, that will be stated below, can be proved by an adaptation 
of the arguments of section 4 of reference [7] ; the only non trivial change is the use of a representation formula 
for a classical solution w of the Neumann problem 

O /* 

(2.31) A a w = 0, in ^ = / on dCl, / aw = 0 

dv Jd n 

(where / is such that a f = 0) which extends the usual layer potential representation for harmonic functions. 
We prove here this formula. 


Lemma 2.9. Let w be the solution to (12.311) . Then 

1 ( 


(2.32) 


Ky) = 


2 tt a{y) 


'an 


a(cr) G a (a,y) f(a) da 


where G a is the Green’s function defined in 

Proof. We first remark that, by the boundary condition in (11.101) and by the normalization condition for w in 
(12.311) above, the following holds: 


/an 


a(a) w(a) d v G a (a, y ) da = 2ira(y)w(y) 


Then, we can compute 

a{a) G a (a,y) f(a) da = / a(a) G a (a,y) d„w(a) da - / a(a) w(a) d v G a (a, y) da + 2na(y)w(y) 
! dan dan 


/an 


= / div 
dn 


G a (x, y) a(x)V x w(x) - w(x) a(x)V x G a (x , y) dx + 2na(y)w(y) 


in L 


G a {x, y) div(a(a;)V x w(a:)) - w(x) div(a(cc)V x G a (:r, y)) dx + 2ira(y)w(y) 

= 2Tra(y) w(y) 

Hence, formula (12.321) follows. 


□ 


As discussed above, one can now reproduce all the estimates proved in [7] (in particular, those in Lemma 4.2 of 
[7]) which lead to the following result 

Proposition 2.10. Let u\ n € H 1 ^), X n -A 0 + be a sequence of solutions to (11.71) given by theorem \l.l[ Then, 
there exists a subsequence, also denoted by u\ n , a regular finite measure m on dfl and a finite set of points 
C dLl, N > 1 such that 

X n | sinh(u^ n ) | —> m 

on dfl in the sense of measures and the points x^ l \ i = 1 are exactly the points at which m has point 

masses. The same points also represent the blow up points for the sequence 

o 1 /' 

U\ n = U\ n - 7 - / au\ n 

Jan a dan 

in the sense that 

= {a: £ Cl : 3 x n —» x, x n £ Cl, with (a-’ n )l —t cxd} 

3. Blowing up solutions 

We now show that there are solutions to problem which concentrate at isolated critical points and 

£ 2 , say, of a constrained on dn as A —>• 0. By the previous discussion, it corresponds to build up solutions to 
problem that concentrate positively and negatively along two geodesics of the boundary of T> which are 
nothing but the Q —orbits of £1 and £2, respectively. 

This section is organized as follows. In subsection 13.11 we write an approximate solution for problem (11.71) . In 
subsection 13.21 we study an associate linear problem and in subsection 13.31 we reduce our non linear problem to 
a finite dimensional one. In subsection 13.41 we study the reduced problem and we prove Theorem 11.31 


1 By an abuse of notation we write here a(o) (and similarly G a ((T,y), etc.) instead of a(x)\. 
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3.1. The approximate solution. To define an approximate solution for the problem ED, a key ingredient 
is given by the solutions of the following problem 

Av = 0 in R 2 


(3.1) 


dv 

dv 


= e on 




where denotes the upper half plane {{x\,X 2 ) : X 2 > 0} and v is the unit exterior normal to 
The solutions of ED are given by 


+ ■ 


W t ^{x !,X 2 ) = log 


2 /r 


(xi - t) 2 + ( x 2 + /i) 2 


where tsK and /i > 0 are parameters. 

Let us provide an approximation for the solution of our problem. Let 

2 /ij 


f{x) = log- 


j = 1,2 £ dCl, n : j > 0. 


I x-Zj -A/ij^| 2 ’ 

In order to satisfy the equation A a u = 0, we need an additional term Hf defined as follows: let Hf(x) to be 
the unique solution of 


—A a Hf = V log a ■ Vu 


Cl 


OH* du A x l 

+ Ae j - A-j:- 

dv dv J a 


“■» on dil 


an a Jon 


id n 


Hf dx = - 


IdSl 


uj dx. 


Now we set 

U\{x) = [u\(x) + h£(x)] - [u%(x) + H$(x)] 
and look for a solution of ED in the form 

u\{x) = U\(x) + $\(x) 

The higher-order term will satisfy some suitable orthogonality conditions (see (13.121) below). 

Let G a (x,y) be defined as in (11.101) and H a (x,y) be the regular part of G a (x,y) defined as in (11.111) . In 
the following we will write simply G(x,y) instead of G\(x,y) and H(x,y) instead of H\(x,y). 

It is immediate to see that H a (x,y) solves the following problem 

x-y 


(A a ) x H a {x,y) = 2V log a • 


\x-y I s 


n 


ddi a , A 2tt 2 (x-y)- v(x) 

-—(x,y) = - 7 - r -a(y) + -,-^- on dCl 

dvx Joa a ( x l F-2/r 

The function Hf can be estimated in terms of H a {x,y). Indeed, the following result holds. 
Lemma 3.1. For any a £ (0,1) 


Hj{x) = H a (x,£j) - log2 Hj + 0( A“) 


uniformly in Cl. 


Proof. The boundary condition satisfied by H A is 

dH j du j , u x A /' 

-JT- = ~^T + Xej - J - / 

dv dv J dQ a J dn 


ae u i 


— 2 A fij 


1 - v(^) ■ v(x) 


(x-£j) -Hx) 


fj \yjv{f,j)\ 2 \x A/ijZ/(^j)| 2 fan a Jdn 


ae< 






















As A -a 0 we get: 

A 
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ae u i 


2 /ij 


j(£j + *hv) 


/an 


X Jon a(x) I* ^ 0 - ViKCi)! 2 2 |y - u(0)\ 2 




= 2 a&) 


an -L |y — v(fi)\ 2 

X^-i 


+ o(\\Va[^)\^ 


aa ~^ I y - K°)l 2 

XHj / 


1 


/a-V " 1 1 + t2 


dt 


2a(Si) \ J-I TTfi d> -° 

a (Cj) (2 tt + O ^arctan(A/ij) _1 — ^)) + O(Xfij) 

2n a(£j) + 0(arctan(A/Lij)) + 0(\pj ) 

27t a(£j) + O(Xnj). 


Let us consider the difference 


Z\(x) = H H X ) ~ H a{x,ij) + log 2 jj,j. 


Hence z\ solves the following problem 


-A a z^ = -A a Hf + A a H a (x,£j) 


dz{ 3 _ dHf dH a 


namely 


dv dv dv 

-A a zi j = Vlog a{x )V [log - log 




dz\ 


= 0(Xnj) 


dv v 3 

As done in Lemma 3.1 of (5] it follows that for any p > 1 

,, dz 


g^\\LP(dU) < C\P. 


Moreover, again as in Lemma 3.1 of [5], we get 


log 


x ~ Cjl 


- log 


1 


Now 


Cy-A/V'(Cy) 

|/i| < C\ 2 ( log 


... = h + I 2 


while for p £ ( 1 , 2 ) 

In conclusion, for any p £ (1, 2) 


L p ( Q) dBx/j.^ (^j )CO dn \B\^^(^j) 

l\p 

A 

\h\ < CA P . 


log 


x ~i> I 2 


- log 




< cA. 


Lp(n) 


By L p — theory 


\ z x Iliv 1 + s .p(n) < C 



dzf 3 


dv 


+ || A q 2 a j ||ip(n) <CX p 


LP(d O) 


for any 0 < s < A By Morrey’s embedding we obtain 

Il4 j Ilc7(n) < cAp 

for any 0 < 7 < | + A This proves the result with a = A 


□ 
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Moreover the function H a (x,y) can be expanded in terms of H(x,y). The following expansion is proved in 
Lemma 2.1 of m- 

Lemma 3.2. Let H ay {x) = H a (x,y) for any y £ f2. Then y —> H a y is a continuous map from 52 into C°' r, (fl), 
for any 7 6 (0,1). It follows 

H a {x,y) = H(x,y) + V log a(y) ■ V(|x - y\ 2 log \x - y\) + H{x, y) 

where y — > %(-,y) is a continuous map from Q into C 1,7 (52) for all 7 £ (0,1)- Furthermore, the function 
(x, y) —> 'H{x , y) £ C 1 (fl x ft), in particular x —>• H a (x , x) £ 


We consider now the following change of variables 


n 


x = Xy, y £ 52 A = —, v{y) = u(Xy). 


Then u is a solution to problem ED if and only if v solves the problem 

in 52 A 

(3.2) 


dv 

dv 


= 2 A 2 sinh v on <952 A , 


where a\(y) := a(Xy). In the expanded domain 52 A , U\(x) becomes 


V(y) = 

3 =1 

where F = X and i/ = i/(£5). Therefore, 


log 


n \ , zj\, 


I y -&-h 


772 ~ 2 log A +i?< (Ay) 




(3.3) 


t>(z/) = V{y) + (!>{y), y £ 52 A . 


will be a solution of (13.21) provided <j> solves 


where we set 

(3.4) 

(3.5) 

and 

(3.6) 


A qa ^ = 0 in 52 a 


^ — W(p = TZ + Af(<p) on dfl\ 
ov 


W(y) := 2A 2 cosliF(y) 


n{y) := - 


ay 9 

—-2A“ sinh V 

ov 


(y) 


Af{4>) = 2A 2 [sinh(I4 + </>) — sinh V — (cosh V)<f >]. 

First of all, we prove that V is a good approximation for a solution to (13.21) provided the parameters yi and /r 2 
are suitably choosen. 

Lemma 3.3. Assume 

(3.7) log2yi = H a (£i,£i) - G a (6,6) and log2y 2 = #0(6,62) - G 0 (6,6)- 

Then, for any a £ (0,1), f/iere exists a positive constant C independent of X such that, for any y £ 52 A , 

O 

1 


(3.8) 

and 


\n{y)\<x a Y, 


- l + \ v-% 


Vy £ 52 A , 


W(y) = J2 1 —r^U + ^(y)), with |0 A (y)|<CA“ + f7A^|y-e'|. 

7=1 (y 67 7=1 


(3.9) 
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Proof. By Lemma l3Jl we deduce that if \y — £)| < j 

~ ( log+(setti,,g *'■=*-&) 

= (H a (Xz + £1, £1) - log2/ii) - 

- (log2 k + log + H a (Xz + {„&) - log2 k) + 0(A“) 

= tfa(a,a) - log2^i - G a (a,6) +o(A a ) + o (\\ z \) 

^-v-' 

=0 because of m 
= 0(A“) + 0(A|i,-£i|) 


and in a similar way 

H^(Xy)~ (log- 7T~~~ - TTTTpT + = 0(X a ) + 0(X\y - ^D- 

V * \y - & - J 

Therefore, the proof follows exactly as in Lemma 3 of [6j. 


□ 


3.2. A linear problem. The key ingredient in this section is the linearization of problem m around the 
solution wo t fj., namely the problem 


A()(> = 0 


(3.10) 

In [6j it has been proved the following result. 


dcf 

dv 


2/i 


on 


14- 


Lemma 3.4. Any bounded solution of (13.10p is a linear combination of the functions 

X2+ y 


and 


Zo(x) = 1 — 2/i 
Zi(x) = —2 


xf + (x 2 + /r ) 2 
xi 


xf + (x 2 + A 4 ) 2 " 


Now, let us assume that the points £1,^2 G dfl are uniformly separated, namely |£i — > d for some d > 0 

which does not depend on A. We have to redefine zo and z\ in a neighbourhood of £1 and £2 in a suitable way. 
So, let Fj : B p (fj) —> N 0 a diffeomorphism, where p > 0 is fixed and N 0 is an open neighborhood of 0 e f 2 
such that 

Fj (ft n B p (£j)) =l(nlVo, Fj (dft n B p {^)) = n No 

and such that Fj preserves area. 

For y € fix H B p / X (f'j) we define 

Fj'i.y) = jFj(Xy) and Z^y) = Zjj(F^(y)), j = 1,2, i = 0,1 
where z l;l denotes the function Zi with parameter a ij, namely: 


z 0j — 1 2a Lj 


X2 + A lj 


xf + (x 2 + A lj) 


2 > 


Zij — 2 


Xl 


xf + (x 2 + Ab) 2 


Let x : K —)• K be a non-negative smooth function with x(r) = 1 f° r x < Rq and x(r) = 0 for r > Rq + 1, 
0<x<l (with Rq a large positive constant). Then, we set 


Xj{y) ■■= x(\Fj X (y)\), j = 1,2 and \{v) ■= Xi{v) + X 2 (y) 


We will assume that A is small enough to satisfy 

\F x (y)\ >R 0 + 1, 


VyeflxndBp/xit'j) 
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Hence, the products XjZij can be defined in the whole domain £l\ by continuation by zero in 
Moreover, by the definition of Z^j we may also assume that, for fixed 0 < b < 1 and suitable chosen <5, 

z 0j (y) > i - A 6 , v y e n x n dB s/x (%) 


We now define: 
(3.11) 


f min (! - A b , Z 0j (y)) if \y - £'-| < f 

z(j/) = ^ 

l 1-A b if|y-,£'|>f for j = l,2. 


We want to solve the following linear problem: given f £ and h £ L°°(dflx)> find (j) £ and 

Cj gR, j = 0,1, 2 such that 


~&aA = f 
dcj) 


(3.12) 


in 


dv 


— W(j> = h + ^ CjXjZij + cqxZ on dtt\ 
i=i 


/ axZ(f> = 0 and / aXjZ\j(f> = 0 for j = 1 , 2 . 


It is necessary to introduce some L°°-weighted norms: if h £ L°°(dfl\) and / £ L°°(Qx), let 

l%)l II r,l _ \f(y)\ 


ft * = sup - 5 - 

^E^ia + ly-^l)" 1 - 

where rr > 0 is a fixed and small number. 

The following result holds. 


and 


** = sup —g- 

^wE- =1 (i + |y-^l)- 2 - 


Proposition 3.5. For any d > 0, there exist Ao > 0 and C > 0 such that for any A £ (0, Ao), for any 
£ 1,62 G dfl with |£i — £ 2 1 > d, for any h £ L°°(dfl\) and f € there is a unique solution (f> £ L°°(Q\) 

and co,ci,C 2 £ K to the problem (13.121) . 

Moreover, 


||0||L~(n x ) < Clog 


A 


+ II/H**) and max{|c 0 |, |ci|, |c 2 |} < C(||/i||* + ll/ll**). 


Proof. We argue as in the proof of Proposition 1 and Proposition 2 of [B] . We only point out 

Aa X t%) = &4>{y) + A </>(;/)» y e o/a. 

a{Xy) 

Moreover, the proof exploits a potential theory argument where Green’s function for the Laplacian is replaced 
by Green’s function G a whose regular part is studied in Lemma 13.21 

□ 

3.3. The non linear problem with constraints. In order to solve our problem we need to split the error 
term in (13.31) as <p(y) = xZ{y) + (j>i{y) where the function Z is defined in (13.111) . r = r(A) is a small parameter 
and satisfies the orthogonal conditions 

/ axZ(f> 1 = 0 and / aXjZ\j(f> 1 = 0 for j = 1 , 2 . 

J £ 1 \ J 

Therefore, the function v in (13.31) reads as 


v (y) = v i(y) + <My)> where V 1 (y) = V{y) + rZ(y) y £ H A . 
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Moreover, v is a solution for (13.21) if and only if (j>\ solves 

( = rVloga A • VZ in 

[ ^-Wi<h=Ki+Mi(<h) on dCl\ 
where (see also ®3D, (ESI) and ESP ) 

(3.13) Wi (y) := 2A 2 cosh V x (y) 


(3.14) 

and 


Ki(y) '■= — 


dV L 9 

—-2A 2 sinh V\ 

av 


(y) 


(3.15) 2Vi(<^i) = 2A 2 [sinh(Vi + </>i) — sinh Vi — cosh(Vi)^>i]. 

It is important to point out that, since Z(y) = 0(1) on all Q\, it follows that V\(y) 

y £ 


V(y) + 0(|t|) for any 


Let us consider first the following auxiliary problem 
-Aa^i = rVlogOA • VZ 


(3.16) 


in fl\ 


d(\> 

dv 


- - yVi(j)i = Ui + Afi(cj)i) + coxZ + ciXiZu + C 2 X 2 Z 12 on dfl\ 


/ a XjZij0-i dx = 0 j = 1,2, 

>. J Cl\ 



axZ(f> 1 dx = 0. 


where Wi, 1Z\ and A f\ are defined in (13.131) . (13.141) and (13.151) respectively. 


Lemma 3.6. Let a £ (0,1), d > 0 and r = 0(A 8 ) with /3 > Then there is Ao > 0 and C > 0 such that for 
any A £ (0, Ao) and for any £±,£2 G dfl with |£i — ^| > d, problem (13.161) has a unique solution (j>\ £ L°°(f2 A ) 
and Co, Ci, C 2 £ R such that 

1101 ||i°°(£ 2 A ) < OA“. 

Furthermore, the function (t,£(,££) —► < Ai( r i £ 2 ) £ L°°(fl\) is C 1 and 

IIAeU^illi'”^) - CXa and ll-^T^i||z,~(n A ) < C\ Pl for some /3i < /3. 

Proof. We argue as in the proof of Lemma 8 of [ 6 ]. The only difference is due to the presence of the R.H.S. 
/ = rVlogOA • Vi? in (13.161) . 

Indeed, first, we point out that 

m(y) = W(y) + 2A 2 sinh (V)tZ + r 2 A 2 cosh(fo + fZ)Z 2 , 

^^ 

:=tB 

where W is defined in (13.41) and |f| < |r|. It is easy to check that ||2?||* < C. Then we write the problem (13.161) 
in terms of the operator A that associates to any </>i £ L°°(fl\) the unique solution given by Proposition 13.51 
with h = T.B 01 +72.1 + and / = rVlogOA • VZ. In terms of A, the problem (13.161) is equivalent to the 

fixed point problem </>i = A{4> 1 ). Therefore, we are going to prove that A is a contraction mapping of the set 

C = {(j)£C{flx) : Mi~ (nj o < A“} . 

From Proposition 13.51 we get 


P(0i)|U~ ( n A) <qiogA| 


MH-B^lll* + ||M(^i)||* + ll^ill* +II/H** 

'--- 

~v 


Arguing as in [6] we get that 


||D||* < C (X a ~ a + A 2/3 A“ +/3 + A 2q ) 
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for some a G (0,1) and a > 0 small so that a — a > a (a is the number in the definition of || • ||*, || • ||** and /? 
is such that r = 0(X /3 )). On the other hand it is easy to check that 


\f(y)\ = 0(Xr\VZ(y)\) 


O | A 1 + V^(l + |y-e'|)- 

3 =1 


for any y G ft \ 


and so 

||/||** = 0(rA a - ff ) = 0( X 0+a -°). 

Then the proof follows exactly as in Lemma 8 of [B]. □ 

Next, we have to choose the parameter r so that the nonlinear problem (13.161) has a solution with Cq = 0. This 
is the result of next lemma whose proof can be carried out exactly as the proof of Lemma 9 in [5]. 

Lemma 3.7. Let d > 0. For any a G (0,1), there exist Ao > 0 and C > 0 such that for X G (0, Ao), and 
an y £i ,£2 G dft with |£i — £2) > d, there exists a unique r with |r| < CX a ~ b / 2 (b is given in (13.111) ). such that 
problem (13.161) admits a unique solution <fi 1 G L°°(fl\), cq = 0 and ci,C2 G K. Moreover 


(3-17) IMk~ (nj o < CA a 

and the function (^(,^2) “t 0x(Ci5^2) C 1 anc ^ 


(3-18) \\D { s U ' 2) <t>i\\ L °°V x) <CX a . 

3.4. The reduced problem and proof of Theorem 11.31 For any (£1, £2) G dfl x dfl with £1 £2, we define 

0(£i,£ 2 ) and Cj(£i,£ 2) for j = 1,2, to be the unique solution to (13.161) with Co = 0 satisfying (13.171) and (13.181) . 
In this section we shall find the points £1 and £2 on the boundary dfl such that ci = C2 = 0. That choice will 
provide a solution to our problem. 


Lemma 3.8. Let a G (5, |) and b G (2(1 — a), 1) (b is given in (13.111) ). It holds true that 


d = tA 


-2 io v ”“ te>+<,<1) 


1,2 


uniformly with respect to (£i,£ 2 ) G dfl x dLl with £1 ^ £ 2 . 

Proof. We multiply equation in (13.16|) by a\XjZij, j = 1,2 and we integrate in y. We take into account that 
V\ = V + tZ and r is choosen so that cq = 0 (see Lemma [XTli . Therefore we get 



r (VaxVZ) XjZijdy 



a\A ax (xjZij)<f>idy 



a\(t>idv (XjZij) dy 



axXjZijdvfady 
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and so 
(3.19) 


2 f 

£ Ci / a(Xy)xiXj z u z ijdy 
,- = i -/aru 


Jo 


= -/ a\A ax (xj z ij)<l>idy - t (Va x VZ) xjZijdy 

+ / axfadv {Xj z ij) dy- a\Xj z ij (WWi + "Ki + M(^i)) dy 
Jdfl\ J d£l\ 

= -/ axA ax {Xj z ij)<l>idy- r(VaxVZ)xjZ lj dy 

J £} x J 12^ 

'-V-' '-V-" 

II 12 

+ [ ax^i [d v (XjZij) ~ Wi (XjZij)] dy 

J dQ\ 


'dClx 


dV 

dv 


— A 2 sinh(V) 


a\Xj z ijdy + r - WZ^j axXjZijdy 


U 


— 2A 2 f axXjZij [sinh(V + tZ) — sinh(F) — rcosh(IA)Z] dy 
JdQx 


Ie 


IdClx 


■Ni(<l>i)axXj z ijdy- 


h 


Now let us estimate each integral TVs of (13.191) . It is immediate to check that 


Indeed if i = j 


We remark that 


27T 

Iq = —a(£i) + o(l) if i = j and Jo = 0 if i ^ j. 


[ a{\y)xiXi z u z ijdy = a(&) f 4Xl 2 <fa+ + o(l) = 2lra ^ + 0 (l) 

ldn x + Hi 


Z 1 j = O 


1 + r 


and VZij = O 


1 + r* 


as r —>• oo. 


We have 


a\A ax (xjZij) = Axi • Zij + 2VxjVZij + x/AZy + VaxVx/Zij + XjVa A VZij 


= O 


A 2 


1 + r / V 1 + r 


A 


+ 0 I , + Xj^-Zij + O —— + O 


A 2 


1 + r/ \ 1 + r 


A 


However, 

and hence 

hence 


A y Zij = A xZl + 0(A|t||V 2 2i|) + 0(A|Vzi|) 


AZ y = O 


1 + r z 


O 


1 + r 


/i| < CAIHI™,) < CA 1+ “ = 0 (A 1+ “-5). 



























18 


CARLO D. PAGANI, DARIO PIEROTTI, ANGELA PISTOIA, AND GIUSI VAIRA 


Moreover, 


< 


1 dClx 


axfadv (Xj) z ijdy 


'dflx 


a\4>iXj {d^Zxj - Wi Zij) dy 


< C\\\(j>i |U~ ( n A ) log \ + C\ a ||0i |Uc. (n ,) log \ 


< C ( A i+ “ log - + A^ Q log — ) = o(A i+a_5 ) 


because Vxj = 0(A) and (see also (3.39) in [B]) 


Xj [dvZu - mzu) = O 


A° 


i + |y-£l 


and since b > 2(1 — a). 
By (13.811 we get that 


'll = (z^ 1 (Ay) - zl 2 (Xy))e H ^ Xy ^ ) - H ^^'> + 0( A 2 ) 


and hence by making the change of variable x = F x (y) and by observing that (F x ) 1 (a:) = x + Cj + 0(A|x|) we 
get 


h = 2 fij J z^{\x + ^j)a{\x-\-ij)e Xx+ °^^ 2) d xl -^ 

/■£ 

= 2 iij d Xl z x (Ax + £j)a(\x + ^)e Ax+0(A2|x|2} 


a + M? 




—2nj\ / d Xl a(^j + Xx + 0(A 2 |ai|))z ^(Xx + £j)e Ax+ °l A ^ 


xi + tf 


Sire, i , on2uu„n„ , t .\„\x+0(\ 2 \x\ 2 ) 


and so 


— 2y,jX J d Xl z x 3 {!;j + Xx + 0(X \x\))a(Xx + t;j)e 

£_ 

- 2/ijA f z x (£j + Az + 0(A 2 |a;|))a(Aa; + t;j)d Xl e Xx+ °( x 1*1 


\h\ < A 1+a log - = o(A 1+a "5). 
A 




xl+tf 


As done in the estimates proved in p. 211, |5 we get that 

Xj {d v Z - WZ) = (zl'iXy) - z b(\ y ))e H *Vv&)-H a (t i &) + 0 (X 2 ) 
and so, by making computations as before we get that 

|4| < OA 1+Q r log i = o(A 1+ “-^). 

A 

Moreover, since by mean value theorem 

XjX 2 [sinh(V + tZ) — sinh(lA) — t cos\\{V)Z] = XjT 2 sinh(V + tZ)Z 2 


by making again the same computations as before we get we also have 

|4| <OAt 2 logi = o(A 1+ “-5). 

A 

Finally, we have 

\h\ < C||0|lioo (Ox) log < CA 2 “ log i. 

Therefore, it remains to estimate the leading term 4 of the R.H.S. of (13.191) . We observe that in the function Z 
is constant except in the regions | y — £' | < yjX ~ b / 2 , j = 1, 2. For sake of simplicity, we can also assume that the 
boundary of 12 in a neighbourhood of the point can be described as a graph of a smooth function ipj defined 
in a neighbourhood of 0 such that ipj{ 0) = ipb (0) = 0 so that 

Fj(si, s 2 ) = (si,s 2 - ^j(si)) and F~ 1 (t 1 ,t 2 ) = (4,4 + 
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Let us remind that F^(y) = Fi ^ v ^ , then 

h = [ t (\7a x {y)VZ(y))xj(y)Zij(y)dy 


= r f Va(Ay)V (z oj (Ffiy))) x (|i?(s/)|) Zij (Ffiy)) dy 


= T 


da\ 

dyi 


dz 0j d(Ff )i dz 0j d(Ff ) 2 


da\ 

dyi 


= tX I 


dxi dyi 8 x 2 dyi 
xx(\F}{y)\)zij{F${y))d y 


Bzoj «(*?) 1 a Zoj «(F/)j 


SoR <9?/2 cAc 2 9?/2 


XX (| F /(y)|) zij (Ffiy)) dy 

t T^\ t \ • F _1 (Arc) x 

(we set rj[y) = i.e. y = —^—) 


= rA 


l (— 

'B(o,p)n*2 l dyi 


(\XI,\X 2 + <Pj(Xxi)) 


dzoj / s. dzoj . / , 

—W- — 


x x(N)*u 


+tA 


/B(0,p)nKl 




{^(Axi, Ax 2 + ^(Aa; 1 ))^^-^)} x x(kl)^ij ( x ) dx 
f d ~ 0j (x)z lj (x) dx + o(l)'] = tA f^-(0) [ (~8p x i ( x 2 + y-j) 


dx 


dyi 


3 { X \ + { X ^ + H y)\ 


dx + o( 1 ) 


= rA (- 7 rVgna(Zj) + o(l)) 


We point out that this is the lower order term of the R.H.S. of (13.191) . because its rate is of order A 1+a b//2 
because of the choice of r. 


The claim follows collecting all the previous estimates. 


Finally, we conclude the proof of Theorem 1 1.31 

Proof. [Theorems 11.31 completed ] Let and £ 2 be two different C 1 —stable critical points of a restricted on <912, 
namely the local Brouwer degree 

deg (Vgna, B(&,p) n <912, 0) ^ 0 for i = 1, 2 

provided p is small enough. Then by the product property we deduce 

deg((Vana, Vona), (B(£i,p) x B(£ 2 ,p)) n (<912 x <912), 0) ^ 0, 

which implies together with Lemma (13.81) that if A is small enough there exists (£i ,£ 2 ) which approaches (£i,£ 2 ) 
as A goes to zero such that ci = c 2 = 0. Therefore v\{y) = V(y) + cf(y) y G 12^, turns out to be a solution to 
the problem (13.21) . It is clear that the scaled function u\(x) = v\(x/X), x G 12 is the solution to problem (11.71) 
which satisfies (11.121) and concentrates at the points and £ 2 as A goes to zero. □ 
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